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We derive a set of equations describing the linear response of the convergence properties of a geodesic con¬ 
gruence to arbitrary geometry variations. It is a combination of equations describing the deviations from the 
standard Raychaudhuri-type equations due to the geodesic shifts and an equation describing the geodesic shifts 
due to the geometry variations. In this framework, the geometry variations, which can be chosen arbitrarily, 
serve as probes to investigate the gravitational contraction processes from various angles. 

We apply the obtained framework to the case of conformal geometry variations, characterized by an arbitrary 
function fix), and see that the formulas get simplified to a great extent. We investigate the response of the 
convergence properties of geodesics in the latest phase of gravitational contractions by restricting the class of 
conformal geometry variations to the one satisfying the strong energy condition. We then find out that in the 
final stage, / and D ■ Df control the overall contraction behavior and that the contraction rate gets larger when 
/ is negative and |/| is so large as to overwhelm \D- Df\. (Here D - D is the Laplacian operator on the spatial 
hypersurfaces orthogonal to the geodesic congruence in concern.) 

To get more concrete insights, we also apply the framework to the time-reversed Friedmann-Robertson- 
Walker model as the simplest case of the singularity formations. 

PACS numbers: 02.40.-k, 04.20.-q, 04.20.Dw 


I. INTRODUCTION 

It has been known that the spacetime manifolds with high 
symmetries often contain singularities. Typical examples are 
the initial big-bang singularity in the spatially homogeneous 
and isotropic universe models and the curvature singularities 
in the black-hole solutions. The singularity theorems [1^] 
have shown that, however, the spacetime singularities do not 
result from high symmetries, but are quite general features of 
spacetimes satisfying reasonable physical conditions, such as 
the energy condition and no closed timelike curves. 

Though the singularity theorems show the generality of sin¬ 
gularities in physically reasonable spacetimes, it is also true 
that their statements are too universal and general to get de¬ 
tailed information on gravitational collapses. 

Among these processes needed to be clarified, the black- 
hole formations are especially important ones. Though the 
spacetime structures after black-hole formations are quite 
well-understood by the uniqueness theorem for the Kerr so¬ 
lution and its related theorems [2, 3], little is clarified about 
the black-hole formation processes themselves. Indeed, the 
cosmic censorship hypothesis [5, 6] and the hoop conjec¬ 
ture [7], which are the two central conjectures for the black- 
hole formations, have not been satisfactorily proved so far 
even though there have been no physically reasonable model 
found manifestly contradicting with these conjectures [8]. 
One of the reasons for this situation might be the fact that 
there is no established framework for analytically describing 
the black-hole formation processes. 
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Here let us pay attention to the Raychaudhuri equation [2, 
3, 9-11], which describes the focusing property of a given 
geodesic congruence. The equation indicates that, once the 
expansion 6 along a geodesic in the congruence gets negative, 
it approaches to within some finite proper-time along the 
geodesic provided that the strong energy condition is satisfied 
(see Sec.VI). This phenomenon signals the occurrence of the 
conjugate point (the focal point) in the future which in turn 
implies the singularity (in the sense of the timelike geodesic 
incompleteness) in the future within a finite proper-time pro¬ 
vided some other conditions are also satisfied [2-4], 

Considering the above fact, we here choose the strategy to 
pay attention to the convergence properties of the geodesic 
congruence as the starting point for the analytical description 
of the gravitational contractions. More specifically, we here 
aim at constructing a theoretical framework describing how 
the convergence properties of a given time-like geodesic con¬ 
gruence are influenced by the slight variation in geometrical 
properties around the geodesic congruence. In other words, 
we study the linear response of the convergence properties of 
the geodesic congruence to the arbitrary geometry variations. 
Here the origins of the geometry variations are not specified 
and they can either be real physical processes or virtual dis¬ 
placements. Since the geometry variations can be arbitrarily 
chosen by hand, then, we might be able to use them as probes 
to investigate the gravitational contraction processes from var¬ 
ious angles. 

We now show the plan for constructing the framework. 
The outline of the physical process in question is as follows: 
When the geometry is varied, geodesics are also shifted ac¬ 
cordingly, which in turn changes the convergence properties 
of the geodesics. The construction of the framework, thus, 
naturally consists of two steps. As the first step, we establish 
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a key equation relating the geometry variation to the geodesic 
shift. As the second step, then, we investigate the variations in 
the convergence properties of the geodesic congruence caused 
by these geodesic shifts. Combining these two steps, we can 
construct a framework which describes the changes in the con¬ 
vergence properties of the geodesic congruence due to the ge¬ 
ometry variations. 

One key point of the above plan resides in the first step. 
Since the geodesic is the integral curve of a tangent vector 
field, it is a global object, which clearly causes some difficulty 
in pursuing our plan. We shall tackle this problem by intro¬ 
ducing a suitable vector describing the geodesic shift. Fur¬ 
thermore it turns out that using solely the component of this 
geodesic shift vector orthogonal to the geodesic is essential, 
making the expressions much simpler as well as giving much 
more transparent interpretations of the equations. In this way, 
we shall find out the key equation. 

Another key point is to introduce a 1 -parameter family of 
geometries on a fixed manifold ^, by means of which we can 
mathematically handle the geometry variations and derive a 
set of equations describing their influence on the convergence 
properties of a given geodesic congruence. 

After establishing general formulas, we shall then pay spe¬ 
cial attention to the case of the conformal geometry variations. 
This case not only makes all the expressions much simpler, but 
also is quite important for several applications. We shall ana¬ 
lyze the influence on the convergence properties of geodesics 
in the final phase of gravitational contractions by restricting 
the class of conformal geometry variations to the one satisfy¬ 
ing the strong energy condition. We then find out that, when 
f{x) is the function characterizing the conformal geometry 
variation, / and D ■ Df control the overall contraction be¬ 
havior in the final stage (where D ■ D is the Laplacian op¬ 
erator induced on the spatial hypersurfaces orthogonal to the 
timelike geodesic congruence in question), and that the con¬ 
traction rate gets larger when / is negative and |/| is so large 
as to overwhelm \D ■ Df\. 

This paper is organized as follows. In Sec.II, we shall re¬ 
view the established results regarding the deformation prop¬ 
erties of timelike geodesic congruences needed in the paper. 
The basic quantities of our analysis - the expansion 0, the 
shear Gab and the twist cOat- and the Raychaudhuri-type equa¬ 
tions shall be introduced there. This section is also served 
to fix the notations and the definitions of the terms adopted 
throughout the paper. The main parts of the paper are Sec.III 
and Sec.IV. In Sec.III, we shall derive a key equation de¬ 
scribing the geodesic shift caused by a geometry variation. In 
Sec.IV, then, we shall derive a system of equations describing 
changes in the deformation properties of the geodesic congru¬ 
ences due to geometry variations. In Sec.V, we shall focus 
on the case of the conformal variations and shall see all the 
formulas reduce to much simpler counterparts in this case. 
We shall then apply our framework to one concrete, explic¬ 
itly analyzable model; the conformal variations applied to the 
time-reversed Friedmann-Robertson-Walker model. Based on 
the results of the previous section, we shall analyze in Sec.VI 
the changes in the contraction properties due to the conformal 
geometry variations in some detail. Section VII is devoted for 


the summary and several discussions. 


II. BASIC FORMULAS 
A. Timelike geodesic congruences 

Let {^,g) be a smooth n-dimensional spacetime manifold 
with signature (—h H-h). 

Let Vfl be the standard covariant derivative on sat¬ 

isfying the metricity condition (Vagbc = 0 ) and the torsion- 
free condition ((V^V^ — ^b^a)f = 0 for any smooth scalar 
function /). 

From now on, any timelike geodesic y is assumed to be 
affine-parametrized by the proper-time T along it. Thus the 
tangent vector of 7 , = {d/dxY, satisfies [ 12 ] 

with = . ( 1 ) 

Let us take a timelike geodesic congruence over an open 
set Q. in it means that we consider some family ‘rf of 
timelike geodesics among which there is one and only one 
geodesic passing through any given point p in Q.. 

We note that, for a given Q., one can choose continuously 
infinite number of timelike geodesic congruences over Q. as 
is easily seen for the case of the Minkowski spacetime for 
{^,g). 

For a fixed timelike geodesic congruence let 

{ 7 s('^)}sG(ii,s 2 ) ^ smooth 1 -parameter family of geodesics 

in 'if’: Then {7s(T)}sg(^j forms a smooth 2-dimensional 
submanifold embedded in with {t,s} providing a system of 
smooth coordinates on it [13]. By suitable parameterizations 
w.r.t. (with respect to) s and suitable synchronizations w.r.t. T 
(i.e. choosing suitably the T = 0 point on each member of the 
family), one can assume that the tangent vector ^“(t) of 7 j(t) 
and the deviation vector rjs‘‘{x), defined by rjs“{x) = {d/dsY, 
are orthogonal to each other for a fixed s. One can thus as¬ 
sume that geodesics in the 1 -parameter family { 7 s(T)}sg(^i s^) 
are so synchronized that all the T = const curves drawn within 
{ 7 s('^)}sG(ii,s 2 ) orthogonal to the geodesics themselves. 

The deviation vector 77 ^( 1 ) along a fixed geodesic 7 is mea¬ 
suring the deviation of the nearby geodesics among the family 
from being parallel to 7 . 

Here we note that, for a fixed geodesic in the family, the 
pseudo-norm of 77 ^, 77 ^ 77 ^, is not constant along the geodesic 
in general, contrary to the case of (Eq.(l)). This fact is 
important since it allows us to describe the convergence prop¬ 
erties of the geodesic congruence in terms of 77 ^ (see the ar¬ 
guments after Eq.(6)). 


B. Expansion B, shear cUfc twist cOah 

Let 7 be some time like geodesic and be its tangent vector 

as before. Then the symmetric covariant tensor hab defined by 

^ab = Sab ~i~ ^a^b i (2) 
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satisfies hab^^ = 0 so that it is regarded as a spatial tensor in 
the following sense. 

Let Tp^ be the tangent space at a point p in 7 . Let Tp^^ 
be the (n — 1)-dimensional subspace of Tp^ consisting of all 
the vectors orthogonal to Then Tp^^ represents all the 
vectors in the spatial directions for an observer moving along 
7 , so that vectors and tensors constructed by Tp^^ are re¬ 
garded as spatial objects for the observer. In particular, when 
hab in Eq.(2) is regarded as restricted on Tp^^, the former 
plays the role of the spatial, positive definite induced metric 
on the latter. 

Furthermore, the (1,1)-tensor hj’, obtained by raising one 
index of hab by satisfies ha^hb'^ = 5a‘^, playing the role of 
the projection operator 

h : TpJK —^ TpJK^ . (3) 

Now, for a given 1-parameter family {7i}iG(sj^i2) ^ time¬ 

like geodesic congruence we can estimate 
the rate of change in the deviation vector 77 “ along a geodesic 
7 (= 7s=o) in {7 s}sg(^i,v 2)- One can easily get [3] 

, (4) 

ax 

where we have introduced 

Bab = '^b^a ■ ( 5 ) 

To get Eq.(4), we have used [ 1 ^, 77 ]“ = V^77^ — = 0, 

which follows due to the fact that and 77 ^ are vectors as¬ 
sociated with the coordinate functions x and s, respectively. 
Here we should note the positions of indices a and b in the 
definition of Bab, describing ^ 77 “ as a linear transformation 
in the conventional form as the R.H.S. (right-hand side) of 
Eq.(4). We also note that, even though we have chosen the 1- 

parameter family { 7 j}jg(^j ^ 2 ) for introducing the lat¬ 
ter describes the deformation properties of as a whole inde¬ 
pendent of the choice of the 1 -parameter family { 7 s}^G(s(,i 2 )- 

It is important that Bab is a spatial tensor, satisfying 

Bab^‘’ = 0 , ^‘‘Bab = Q , 

as is shown with the help of Eq.(l). 

The tensor Bab plays the central role in our analysis below. 
It is convenient to decompose Bab into the following form; 

1 

Bab — 7 ^Itab T XJab “b (X)ab • (6) 

n — 1 

Here the first term on the R.H.S. is the trace-part of Bab 
with 0 = Ba^', the second and the third terms are both trace- 
free with the former and the latter being symmetric and anti¬ 
symmetric in the indices, respectively. Accordingly it follows 
that Oba = (^ab with = 0 and (Oba = -(Oab- 

It is conventional to call 0, Gab and (Oab the expansion, the 
shear and the twist, respectively. 

Due to Eq.(4), which is in the standard form describing 
deformations in the linear transformations, one can convince 
oneself of the interpretations for 0, Gab and cOab- For a given 


timelike geodesic congruence (synchronized as discussed 
in Sec.II A), the set of the points with the affine-parameter 
X forms an (n — 1)-dimensional section of When com¬ 
paring and then, 0 , Gab and (Oab are interpreted 

as the rate of the dilation of the (n — l)-volume, the shape- 
change and the rotation, respectively, averaged during the 
time-interval At. 

By the Frobenius’s theorem, (Oab = 0 holds iff the timelike 
geodesic congruence is hypersurface orthogonal, i.e. the 
whole of if can be foliated by smooth (n — 1 )-dimensional 
spatial sections orthogonal to every timelike geodesic con¬ 
tained in if [3]. In such cases, Bab reduces to a spatial 
symmetric tensor, coinciding with the extrinsic curvature Kab 
for the smooth spatial sections orthogonal to all timelike 
geodesics in At the same time, (Oab = 0 holds iff the spa¬ 
tial covariant derivative Da is torsion-free (see Eq.(A12) and 
arguments following it in Appendix A). Here Da is the spa¬ 
tial derivative operator on T induced from through the 
projection h (Eq.(3)), defined as 

DaTat = haPh\ ■ ■ ■ h/■ ■ ■ VpTf :/ (7) 

for any tensor Tfff 


C. Notations 

At this point, let us summarize the notations employed be¬ 
low. 

(i) We often adopt the standard notations for the products 

among matrices and vectors when they save indices. 
For instance, second-rank tensors and 5“^ can be 
regarded as matrices. Then 

{AB)\=A‘‘aB‘^b , {Auy=A‘‘bu'’ . 

We also allow the remaining indices to be raised or low¬ 
ered freely as usual tensorial indices. For instance, 

{AB)ab=AacB^b ■ 

(ii) For vectors u“ and v“, we define 

M • V = UaV“ . 

Similarly, for tensors Gab and j5ab, we define (note the 
index-positions) 

(ct • l^')ab = (tCacj^b , 

and 

a-j3 = (a-j3)/ = a«ij3"'’ . 

Note the difference between (a • I5)ab defined here and 
iali)ab = a/^cb defined in (i). 
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(iii) Overlines (underlines) attached to more than two indices 

indicate symmetrization (anti-symmetrization) of these 
indices. For instance, 

(hh= 2+ acbad) , 

(^ab = 2 l^ba^ i j^ab^d = {l^abcd l^cbad) • 

(Anti-)symmetrizations are operated prior to contrac¬ 
tions when they appear simultaneously. For instance, 

g‘^^{rBU = ^g‘‘^{rcdB\-YbdB^) . 

(iv) Any index a is replaced by a dot ( “ ■ ” ) when it is con¬ 

tracted with a geodesic tangent vector For instance, 

Ua.=U,be , V-=V% , r:=Yab^‘^e • 

When no confusions are caused, these dots may often 
be omitted. For instance, 7 ^ can be used as a shorthand 
for 7 “' provided that it is obvious in the context. 

(v) A second rank spatial tensor with the tilde-symbol (“ ~ ”) 

indicates its trace-free part. For instance, when Uab is a 
spatial tensor, 

t^ab — tXab ^ tXc hab ■ 

n — 1 

Furthermore, (a • a)ab is also a spatial tensor. Then, 

{a-a)ab = {a-a)ab - ^—{a-a)hab ■ 

n — 1 


(vi) The projection map Eq.(3) induces maps from tensors to 
their corresponding spatial tensors. These maps are de¬ 
noted by the “underline” symbol attached to the main 
letter representing the mapped object. For instance, 

Ya = haVb , Y^i^ = h/hb'^Ycd ■ 

We also note that, for instance, unambiguously in¬ 
dicates the trace-free part of the spatial projection of 
Rab, noting that “ ~ ” is meaningful only for the spatial 
objects ((v)). 


D. Raychaudhuri-type equations 

Now, to estimate the rate of change of Bab along a geodesic, 
one takes its derivative w.r.t. T. Noting that ^ and 

taking into account Eq.(l), it is then straightforward to get 

-J-Bab =—{BB)ab —Ra-b- , (8) 

dT 

where we adopt the standard definitions for the curvatures [3], 

{'^a'^b-'^b'^a)Uc = Rab/ud for Vm“, 

Rab = Racb , ■ (9) 


We note that -^Bab is a spatial tensor as is seen by Eq.(8). 
In general, -^Aabc - is spatial when is spatial, which is 
easily shown as 

^‘‘^Aabc- = ^{^‘‘Aabc-)=0 . 
dx dx 

By decomposing the both sides of Eq.(8) into the trace-part 
and the trace-free part, and by decomposing the latter into the 
symmetric and anti-symmetric parts, we finally obtain 

de 1 7 

— =- - d - a■ a-Y (o-(o-R.. , 

dx n—\ 

-7- =- - e aab-{(y-Cy)^i, + {(0-(0)ab-Ra.b. , 

clT n — i 

^ = -^ecOab-2{co-a)ab , (10) 

dx n — 1 — 

where Ra-b- is the trace-free part of Ra-b- regarded as a 
second-rank spatial tensor (see (iv) and (v) in Sec.IIC). Eur- 
thermore, it is easy to see that Ra-b- is expressed as 

Ra.b.=Ca.b. - XtKab ; (H) 

n — 2 

where Cabcd is the Weyl tensor. 

The first equation in Eq.(lO) is often called the Raychaud- 
huri equation [9-11], which plays the central role for proving 
the singularity theorems. Eor brevity let us call the three equa¬ 
tions in Eq.(lO) the Raychaudhuri-type equations. 


III. DEVIATIONS OF RAYCHAUDHURI-TYPE 
EQUATIONS DUE TO GEOMETRICAL VARIATIONS 

A. One-parameter family of geometries on .JY 

Let us consider a 1-parameter family of n-dimensional 
spacetime geometries on {{-^,g^^'^)}\eA- Let be 
the covariant derivative compatible with g\i, ■ 

Eor some quantity A (A) regarded as a smooth function of 

O 

A, let A indicate the derivative of A (A) w.r.t. A at A = 0 and 
call the o-derivative of A (A) for brevity. 

Eor notational simplicity, we shall write just A for A(0), 
such as V and gab for and gj^J, respectively. 

We here introduce two quantities in connection with the o- 
derivative for later analysis. 


Y^b=Lb , ■ (12) 

The former represents the geometrical variation while the lat¬ 
ter describes the shift of a geodesic due to the geometrical 
variation. 

By taking the o-derivative on both sides of gabg^‘^ = Ba‘^ 
^a=gab^^, we also get 


lab 


-f‘’ , la = Xa + Ya 


(13) 
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Taking the o-derivative on both sides of = — 1, we get a 
useful formula 

= , (14) 

where Eq.(12) and Eq.(13) have been used. 

Now, according to general properties of covariant deriva¬ 
tives [3], for any vector Wa, the difference between VaWb and 

should be represented as '^‘^ab{^)'^c with '^^ab{^) bC" 
ing some function of X, 

= VaWfo-'^^a&(/l)wc . (15) 

(A.'l (A.'l 

Applying Eq.(15) to the identity V), ’= 0 and following 
the standard procedure, ^‘^ab{^) can be expressed as 

= . (16) 

Here is the covariant derivative compatible to gab so that it 
is independent of X. 

Taking the o-derivative on both size of Eq.(16), then, it fol¬ 
lows that 

(Vfl7M + '^bYad - '^dYab) ■ (17) 

Erom Eq.(17), we get a useful formula 

^‘^ab^‘‘^‘’ = -^7^-U^^7abK‘^^‘’ ■ (18) 

ax 2 

B. Key equation for geodesic shift due to geometrical 
variations 

When a spacetime geometry is varied, geodesics should be 
shifted accordingly. We now derive a key equation describing 
such a geodesic shift caused by geometry variations. 

Consider a geodesic in the spacetime {^,g^^')). Its 
tangent vector, ^^(A), should satisfy 

^*(A)VW<^«(A) = 0 . (19) 

Taking the o-derivative on both sides of Eq.(19), we get 

which along with Eq.(18) yields 

^X‘‘+B\x‘’ = -^7‘^+^^{^'^7bc)e^^ . (20) 

What we are concerned with is, however, the deviation of 
the geodesic from the original geodesic 7 caused by geom¬ 
etry variations (Yah), so that it is more appropriate to use 
the component of x“ orthogonal to rather than x^ itself. 
Indeed, it turns out that the formulas below become much sim¬ 
pler and more transparent in meaning in terms of x±‘‘ [14]. 
On the other hand, the component X\\‘^ which is parallel to 


does not shift 7 , so that it may be regarded as the “gauge- 
freedom” in the present description. In effect, we fix the 
gauge-freedom by Xif = 0 . 

Let us then introduce 

v" = Xx'’=r , ( 21 ) 

which satisfies with the help of Eq.(14) 

v“=z“-^7..r, 

^v = <^~=0. (22) 

at 

Considering Eq.(12) along with Eq.(22), we note 

l" = v“ + ^7..<^" , 

° 1 

^a = Va + 7a + :^7-^a ■ (23) 

Now Eq.(20) can be represented in terms of by means of 
Eq.(22); 

— v‘^+B\v‘’ = -D‘‘Y..-—f , (24) 

dx 2 dx- 

where Da is the spatial derivative operator induced from 
(Eq.(7)). 

It is obvious that Eq.(24), the relation solely among spatial 
quantities, is much more desirable than Eq.(20). 

Equation (24) is our key equation, telling us the linear re¬ 
sponse of the geodesic shift v" to the geometry variation jab- 
Equation (24) is also expressed as 

= , (25) 

2 dt- 

with 

= + . (26) 

dx 

Since (rather than j") contains the pure geometrical infor¬ 
mation on the geodesic shift, it should hold = 0 Yab = 
0, so that should be invertible. Thus one can formally 
solve Eq.(25) as 



which formally expresses the geodesic shift in terms of the 
geometry variation jab- 

IV. CHANGES IN GEODESIC CONVERGENCE 
PROPERTIES DUE TO GEOMETRY VARIATIONS 

A. Deviations of 0, dab and due to geometry variations 

Now taking the o-derivative of 
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we get after some straightforward calculations, 

Bab = (7B)a5 + \y..Bba + \ ^7^1, + DyVa 

+ iBv)a^b + ^a{vB)i, + DbYa + ^a{7B)b ■ (28) 

Here we note that, according to Sec.IIC, all the terms are 
unambiguously defined in shorthand notations. Namely, 

(y^hb = TacB^b ’ = BacV '', {vB)b = v‘^Bcb , 4(7B)fo = 

^aYcd^^B'^b and so forth. 

, O 

Contrary to the case of ^Bab (see Eq.(8)), Bab is not spa¬ 
tial any more as is seen in the second line on the R.H.S. of 
Eq.(28). This result turns out to be a reasonable one, repeat- 

O 

ing the similar argument after Eq.(9) along with 0. 
Similarly we also get 

B% = -g^^ijBU + iy.. + i + DbV‘^ 

+{Bv)%b + i‘'{vB)b+g‘^^{Db7c_+U7_B)b^ ■ (29) 

Here the difference between the first terms in Eq.(28) and 
Eq.(29) should be noted. (Note the signs and the (anti- 
)symmetrization.) 

By contracting both sides of Eq.(28) with g"^ along with 

O 

simple manipulations, we get the formula for 0 as 

0 = ^7--0 + ^^7 +-D-v , (30) 

where y is the shorthand notation for h^^jab, the trace of y^b- 
Equation (30) is one of our basic formulas describing changes 
in geodesic convergence properties caused by geometry vari¬ 
ations. 

O O 

To get similar basic formulas for a ^b ™d (O ab, some more 
considerations are helpful. Eirstly we note that 

O O 

^ ab ^ ba ~ {^ab ^ba^ ~ 0 j 

ttt ab^ ^ ba ~ {^ab T ^ba) ~ 0 j (31) 

o o 

due to the linearity of the o-derivative. Thus, O’ ab and CO ab 
are symmetric and anti-symmetric in the indices, respectively 
just as Gab and cOab- 

Now taking the o-derivative on both sides of Eq.(6), we see 

O O O 1 O 1 ^ 

^ ab~^ ^ ab ab 7 0 - 0 h ab • (32) 

n—\ n—\ 

o 

Considering Eq.(31), then, we can get the formulas for a ab 

o 

and (0 ab by decomposing the R.H.S. of Eq.(32) into the sym¬ 
metric part and the anti-symmetric part, respectively. Taking 


into account Eq.(28), Eq.(30) and Eq.(A13), we finally obtain 

O 1 1 - / \ 

cy ab = :^y..cyab + (^^ab + {y(y)ab + ( 7 ®) ab 

+ ^ -^fab + ^b + 2&((^ v)y 

o 1 

(Oab = -:jy..(Oab- 2^a{(OV)^ - DgVh 

+4(;j-^7fo + ((^7)fo-(®7)fo) , (34) 

where the second term on the R.H.S. of Eq.(33) is the only 

O 

trace-part of a ab, which vanishes for the case of the confor¬ 
mal geometry variation (yab hab)- The last lines of Eq.(33) 
and Eq.(34) also vanish for the case of the conformal geome¬ 
try variation ( 7 “ =0). (The case of conformal geometry vari¬ 
ations shall be analyzed in detail in Sec.V.) 

Let us note that cOab = 0 iff the timelike geodesic congru¬ 
ence in question allows smooth (n — 1 )-dimensional spatial 
sections (see Appendix A for more details). Once (Oab = 0 is 
satisfied, it is satisfied all the way along the geodesic 7 as is 
seen by the third equation in Eq.(lO). When discussing gravi¬ 
tational collapses, thus, it is mostly assumed cOab = 0 , since we 
are usually interested in gravitational collapses started from 
ordinary, mild initial conditions satisfying COab = 0 . 

When considering the geometry variations, it is thus a rea¬ 
sonable assumption that the class of variations we consider 
retains the property COab = 0. As discussed above, it cor¬ 
responds to restricting the variations within collapsing ge¬ 
ometries started from ordinary, mild initial conditions. Erom 
Eq.(34), this condition for the geometry variation yab should 

O 

hs CO ab = 0, yielding 

Dgyb - 4 ^ ® 

With the help of Eq.(27), then, Eq.(35) is understood as the 
condition restricting the geometry variation yab to those that 
preserve the well-behaved nature of the geodesic congruence 
described by COab = 0 . 

B. Deviation of ^Bah due to geometry variations 

Noting that 

RacbH^)=RaJ-2Wa%b‘'{^)+2%b'{X)%/{X) , 
we get 

R«c/ = -2V^<|’‘',^ = -2V,r'^,^(7) , (36) 

where 

'^cabiY) = 2 ^^ ^cYab) i 

r^a&(7) = ^g‘'‘‘ydabi7) ■ 
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Now taking the o-derivative on both sides of Eq.(8) with gab 
replaced by g^^j^, we get 



O 


-BacB'^b - BacB% - 
-Racb‘^l%d-Racb‘^^‘'ld ■ (37) 


With the help of Eq.(23), Eq.(28), Eq.(29) and Eq.(36), we 
can rewrite Eq.(37), yielding 


^^Bab] 


2 c/t’ 


ab 


-BaciDbV'' + ^biBvf) - Bab{D‘^Va + ^aivBf) 


-Racbd{2V^‘‘ + ^^f) 

- {{YB)ac + ]^7..Bc}j B‘^b-Ba‘' (^{rB)bc + ^r-Bbc^ 

+ (b,D*y* + ^*{yB)*\ ■ (38) 


Here (see Sec.IIC (iv)); = 

(AB + ZM)^^, the anti-commutator between A and B, so that 

< B, D*Y*~Y ^*iYB)* r is the anti-commutator between Bab 

and DaYj^ -f ^a{YB)b- 

With some computations, it is straightforward to get 


Combing Eq.(38) with Eq.(39), we finally obtain 


d 

1 c/2 


1 f dY 


= ^-X^Yab-,,\-fZ^B\ --\—y..\Bb, 


2 dx 


1 


ab 


2\dx‘ 


2c/t2 

- (BacDbV‘' + DcVaB%) - {^a{vBB)b + {BBv)a^b) 
-Racbd{2v%'^+^^Y^) 


^DaDbY. 


+ iYcd + Y-hcd)B‘^bB\ + ^ ((7sc + Y-hac)B‘^^b) 

-{YBB)a^b- -^DdYb 

- ({YB)dc + ^Y-Bca^ B^b-Ba {{YB)hc + ]^Y-Bb}j 
+ |b, + \ {b{yB)) a^b-\^a {yBb'^ b ■ (40) 


Here we note that {YBB)a = YbB^cB’^a and ^ 5 ( 76 )^ a = 

BabYcB‘^^, faithfully following the notation rules introduced in 
SecJIC. 

We note that {-^BabY as well as Bab is not a spatial tensor. 

In the case of the conformal geometry variation, Yab = 
2fgab, only the first three lines on the R.H.S. of Eq.(40) re¬ 
main. 

To get explicit expressions for (^O) , {'^^ab) and 
( j^fOabY from Eq.(40), one can follow the similar procedures 
to get Eq.(30), Eq.(33) and Eq.(34) from Eq.(28). However, 
Eq.(40) might be enough for the moment and we shall give 
the explicit expressions for the case of the conformal geome¬ 
try variations below. 


2V,r.,,(7)^^ 

1 c/2 I f d \ „ 1 ^ 

“ 2 " 2 ■ j ■ 

+ {Ycd + Y-hcd)B\B‘^^a + ^ (( 73 c + Y-hac)B‘^'Y) 
-2DaYcB‘'-^-YcDaB‘^b-Yc 
-{YBB)a^b-^DaYb ■ (39) 


V. CONFORMAL VARIATIONS TO THE GEOMETRY 

A. Basic formulas for the conformal variations to the geometry 

Erom now on, we confine ourselves to the cases of the con¬ 
formal variations to the geometry described as 

8iY{x)=^^^^^’‘^8ab{x) , (41) 

where f{x) is an arbitrary smooth function. 

According to the first equation in Eq.(12), we then get 

Yabix) =2f{x)8abix) ■ (42) 


It is worth noting that, in case of the conformal geometry 
variation described by Yab = 2f8ab if is smooth function), 
it follows that jab + 7 ■ hab = 0 and 7 ^ = 0 so that all the lines 
but the first one on the R.H.S. of Eq.(39) vanish and the ex¬ 
pression gets greatly simplified. 


Based on Eq.(42), we further get 

Yab ~ 2fhab ? 7 “ 7^ =2{vi 1 )y , Yab ~ 0 

Ya=2f^a , 7- =-2f , 7fl = 0 , 

Yab^Y-hab=Q ■ (43) 



Then, the key equation Eq.(24) reduces to 


dv" u 

dx 

(44) 

or, in the form of Eq.(25), 


j/7v* = -D 7 , 

(45) 


where is given by Eq.(26), which is invertible as dis¬ 
cussed just after Eq.(26). Thus Eq.(44) or Eq.(45) can be for¬ 
mally solved for as 




dx' 


-1 + B 


Df 


(46) 


Erom Eq.(50), we can extract the formulas for (^0) , 

{j-.OabY- 

Eirst of all, we note that 


dd 

dx 


O 


d 

dX |A=o 




=8 


dx 



Noting that g = —2fg^^ due to Eq.(13) and Eq.(42), the 
first term in the last line reduces to —2/^. With the help of 
Eq.(lO) and Eq.(50), then, we get 


Now, when jab = Xfgab, Eq.(28) reduces to 

O f 

Bab = fBab + B)}jVa + —^ab + {Bv)a^b + ^a{'^B)ij . (47) 
dx 

o o o 

Erom Eq.(47) one can extract equations for 0 , a ab and ab as 
before, which of course coincide with the equations reduced 
directly from Eq.(30), Eq.(33) and Eq.(34) for the case Jab = 
2fgab- In order to reduce Eq.(33) to the corresponding equa¬ 
tion (Eq.(48) below), for instance, we note Eq.(43) along with 

7 - (7 = 2 /a/ = 0 , (ya)^^ = 2 /a^fo, ( 7 ®) = 2 /a^j = 0 and 

so on. We thus get 

O {jf 

9=—fO + {n—1)— +D-V , 

o ab= fOab + D-aVj; + 2 ^^(av)^ , 

0)ab=f(0ab-DaVb-2^a{(0V)^ . (48) 


2/ 

dx J n — 1 



2 

n — 1 


D V 


9 


d^f 

+ {n-l)-^-D-Df + 2fR..-2R.-v 
+2fa a-2fco (o-2a Dv-2(o Dv .(51) 


By the similar argument as in Eq.(31), we see that (7^*) 

and {-^(OabY are symmetric and anti-symmetric in the in¬ 
dices, respectively, due to the linearity of the o-derivative and 
Thus {^(Oab) is manifestly trace-free. Eurthermore 
is also trace-free in the conformal cases, which is 

shown as 




-g‘'^^(yab 

dx 


= 2fg‘^‘’-^a,b = 0 , 


o 

As for a ab^ comparing Eq.(33) and Eq.(48), we see that the 
only trace-part on the R.H.S. of Eq.(33), <yhab, van¬ 

ishes in the conformal cases, so that the latter equation be¬ 
comes completely trace-free. We also note that the last lines 
in Eq.(33) and Eq.(34) vanish in the conformal cases, due to 
7a =0. 

Eurthermore the condition for preserving (Oab = 0 (Eq.(35)) 
reduces to a simple form 

DaVb = 0 , (49) 

which means that the geodesic shift should not form any 
rotation in the (n — 1 )-dimensional section orthogonal to the 
flow-line. A typical situation is that is described as a gradi¬ 
ent of some scalar function as is the case analyzed in the next 
subsection Sec.VB (see Eq.(57) below). 

As for the equation for {-^Bab) , Eq.(40) also simplifies 
drastically for the conformal cases. Indeed, only the first three 
lines on the R.H.S. of Eq.(40) remain, yielding 



d^f df df 

= DaDbf -X—BabX -^Bha 

- [BacDhV'^ + DcVaB\^ - (^^a{vBB)b + {BBv)a^b^ 

-2RacbdV^'^ ■ (50) 


where we have used the fact that -^(Xab is trace-free (as is seen 
byEq.(lO)). 

Thus, by extracting the trace-free part of Eq.(50) and divid¬ 
ing it into the symmetric part and the anti-symmetric part, we 
get 

f d \° df - - 

29 - - 

- - 2U(0(0V)-^ 

-X{Racbd--Rcdgab)v^^‘‘ , (52) 

n 

and 

fd \° df 
20 

+-^^Dgyb + {a, £>* - [cu, D-v-f\^b 

40 

+- 7^a{0}V)b+2.^a{{(y,0}}v)b ■ (53) 

n — 1 “ “ “ - 

We see from Eq.(53) and the third equation in Eq.(48) that 
one can set (Oab = 0 provided that the “twist-less” condition 
Eq.(49) is satisfied for the geometry variations. 
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We note that the arbitrary function f{x) determines the ge¬ 
ometry variation jab (Eq.(42)), which causes the geodesic shift 
(Eq.(46)). Then the two sets of equations (the equations in 
Eq.(48) and the equations Eq.(51)-Eq.(53)) describe the linear 
responses of 0, a and (O and those of -^Oab and ^COab to 
the geometry variation caused by the arbitrary function /. 

Since the function f{x) is arbitrary, then, we can analyze 
in detail the relation between the geometry variations and the 
changes in the focusing properties of the geodesic congruence 
with the help of these sets of equations. 


B. Application to the time-reversed 
Friedmann-Robertson-Walker model 

To get some concrete insights for the present framework, 
we here apply it to the case of the time-reversed Friedmann- 
Robertson-Walker (FRW) model. The FRW model describes 
the expanding universe with an initial singularity, so that its 
time-reversed version can be regarded as a simple spacetime 
model describing a gravitational collapse with a final singular¬ 
ity. (See Appendix B for the basic results of the FRW space- 
time.) 

The metric for the FRW spacetime is given in the co¬ 
moving coordinates as 

ds^ = —dt^ + a^{t) {r)dr^ + r^d'&^+ r^ sin^-i^dip^) , 

(^ = -1,0,1) . (54) 

Let us choose as a geodesic congruence ‘rf the set of stan¬ 
dard geodesics compatible with the co-moving coordinates, 
described by x‘‘{t) = ^(f ro tSo ^o) with tq, tJo and (Po being 
some constants. Then in the present frame of coordinates it 
follows 


^“ = ^(1 000) , ^, = ^(-1 000) , 
/zafo = diag(0 a^F sin^tJ) . 

It is then straightforward to get 

0 1 d 

Bab — ^ab — 'xSabfl ~ ~nab ■, 

1 a 


which yields 


0 = 3- , aab = COab = 0 . (55) 

a 

From Eq.(B6) and Eq.(B8) in Appendix B, we see that 

Ra-b- = - Kb , Ra-b-=0 , R..=-3- . 

a a 

As far as d < 0 (a time-reversal invariant relation), thus, it 
follows that R >0. 

In the context of the generic condition for the singularity 
theorems [2, 3], which is stated as “each timelike geodesic 
contains at least one point at which Ra.b. ^ 0”, the trace-part 
of Ra.fo. (i.e. R..) is non-zero rather than the trace-free part 


{Ra-b. )■ Because of this fact, ^ < 0 is guaranteed even though 
(7at, = 0(seeEq.(10)). 

Then Eq.(44) in this case becomes 


v« + -v« = -D7, (56) 

a 

which can be explicitly solved in the form of Eq.(46) as 

WKx)^f^^dtK (57) 

Jto a{t') 

Here we should note that the operator D“ is t-dependent 
through a{t), which is clearly seen in the co-moving coordi¬ 
nates as 


D« = ^(- 


P-{t)F^{r) ^ K{t)r^ ^ fl^(f)r^sin'0 
Thus Eq.(57) is expressed in a more explicit form 

= , (58) 

a[t) a[t) ^ ' 


where D[\f is with a{t) = 1 so that it is f-independent. 
(These two operators are related by D[\\‘ = ) One 

can then easily check that Eq.(58) is the solution for Eq.(56). 

For simplicity, we shall set v^(fo,x) = 0 below. Then 
Eq.(48) becomes 


O Q f 

0=—3-/-I-3- - a{t)D ■ D\j/{t,x) , 

a dT 

O ab ~ tl[t^]DaD-i^\l/(^tjX) , 

^ab = 0 ■ (59) 


We note that Eq.(59) describes how the focusing of geodesics 
starts to deviate from the one in the exact FRW spacetime: 
The conformal geometry variation induces the variations in 
0 and Gab- In particular Gab obtains a non-zero value even 
though Gab = 0 at the beginning. Once Gab gets non-zero, the 
first and the third terms on the R.H.S. of the second equation 
in Eq.(48) also cause the variation in Gab- On the other hand, 
(Oab = 0 is preserved as far as the condition Eq.(49) is satisfied. 

Furthermore Eq.(51), Eq.(52) and Eq.(53) become 


dt) dt^ 


a at ^ 


d 

a 




-F2dD ■ D\^ — D ■ Df . 

(^) =-^bf + ^'^^b^ 

+M(;) 

^y=o. (60) 

dx j 


Here let us pay attention to the first equations in Eq.(59) 
and Eq.(60) for d < 0, d < 0 and / < 0; 
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fdey dy^\d\df 




D-Df{t',x) 


a{t') 


\dT) ^ 


+ 3 


dx^ 

2\dit)\ 


-6 


a dx 
‘D-Df{t',x 


dt' , 

-}i/i 


i{t') 


dt'-D Df . 


The function f{t,x) can be chosen arbitrarily. The above 
two equations serve as concrete examples for the general cases 
of the conformal variations, which is analyzed in the next sec¬ 
tion. 


VI. INFLUENCE OE CONEORMAL GEOMETRY 
VARIATIONS ON THE FOCUSING PROPERTIES OE 
GEODESIC CONGRUENCES 


We now investigate the response of the gravitational focus¬ 
ing properties of geodesic congruences to conformal geome¬ 
try variations. For this purpose, let us focus on the behavior 
of the expansion 9. 

We first recall how 9 behaves during the latest phase of 
gravitational contractions. We pay attention to the Raychaud- 
huri equation, the first equation in Eq.(lO). As discussed at 
the end of Sec.IIB, one can set cOat = 0 when the timelike 
geodesic congruence is hypersurface orthogonal, which is 
usually assumed. Furthermore R.. > 0 is guaranteed assum¬ 
ing that the strong energy condition is satisfied for the matter 
content. Thus 

d9 1 , 

+-7 0<O. (61) 

at n — I 

which is equivalent to 

d9-^ 1 

->- , 

dt n — I 

or 

0(t)< ^0(To)^' + ^^(t-To)^ 

Thus once 0 becomes negative at some To (i.e. 0(to) < 0), 
then 0 ——oo as the proper-time tends to some finite value x 
satisfying To < T < To + [3]. 

The function f{x) is arbitrary as far as the consistency of 
arguments are retained. In the present context of gravitational 
collapses, it is reasonable to enforce the strong energy condi¬ 
tion for the consistency with the argument just before Eq.(61). 
Indeed, taking into account the energy condition for restrict¬ 
ing the class of functions to be considered is sometimes vital 
for describing the gravitational contractions correctly [8]. 

Let us then pay attention to the strong energy condition, 
translated as R.. >0, and derive a formula for {R..)°. By the 
conformal transformation Eq.(4I), the Ricci tensor is trans¬ 
formed as [3], 

-(-(«-2)VaV*/-Ag„£,A/ + 6)(A^) . 


On the other hand, the tangent vector of a geodesic is trans¬ 
formed as 


y^^‘^ = {l-Xf)y + Xv‘^ + 0{X^) , 

where Eq.(23) and Eq.(43) have been used. Thus we obtain 

{R..)° = -2fR..+D-Df-9^-in-l)^ + 2vR. , 

(62) 

where Eq.(Al 1) has been used and v R. = v'^Rab^ (follow¬ 
ing the notation rules in Sec.II C). 

To retain the strong energy condition, it is reasonable to 
choose the function f{x) which does not greatly depart from 
the condition {R..)° >0. 

Now we pay attention to the first equation in Eq.(48), and 
Eq.(51). We consider the final phase of gravitational contrac¬ 
tions where 0 is negative and 101 At ^ 1 according to the ar¬ 
gument after Eq.(61). Here At is the typical time-scale of 
changes of our concern. Let Atq be the time-scale in which 
0 changes substantially, then the above estimation implies 
Atq ^ At. Next, the behavior of f{x) is quite arbitrary but 
at least Atq <C AT^ ^ At and Atq ^ AX^ ^ At should sat¬ 
isfy in order for our arguments to be meaningful. In particu¬ 
lar, At, the typical time-scale of our concern, should be large 
enough compared to the other scales to investigate the effects 
of geometry variations. Thus we impose 

At0 < (A7> , KKf) < At , (63) 


implying that 

l<l»(lgl, IO-I3/I). 

Eor simplicity, let us assume / < 0 also, which roughly cor¬ 
responds to the conformal geometry variations in favor of the 
geodesic focusing. 

Now we focus on the first equation in Eq.(48). Let us first 
analyze the importance of each term on the R.H.S. of the equa¬ 
tion. We see from Eq.(46) (along with its counter-part for the 
time-reversed LRW model, Eq.(57)) that 

£).v~-AtD-£)/~--^/ . (64) 

AXf 

Then we estimate as 

l-P-v| AtAtq I^I Axe 1^1 AX/ 

|/0| ^ AX/ ’ ATf ’ |£).vrA7>AT ■ 

There is no a priori relation between ATf and AXf. One rea¬ 
sonable situation is, however, that they are of similar mag¬ 
nitudes, e.g. ATf ~ AXf, unless quite special situations are 
considered. 

Thus, along with Eq.(63), we see 

(|/0|, |r>-v|) » igl , (65) 

implying that the term £) • v is as important as the term /0 
and they are much more important than the term ^. 
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Next let us search for the condition for 0 < 0 which implies 
the stronger focusing of geodesics due to the geometry varia¬ 
tions. In later phase of the geodesic congruence contractions, 

O 

the threshold condition 0=0 becomes 


dx 


|e| 

n—\ 


■/= 


1 


n—\ 


D V , 


which along with /(Tq) = fo can be solved as 


( 66 ) 


m = (/o - ■ V , 

A(t,To) = [ |0(T')|<iT' . 

n - 1 Jto 


We note that 


A(t,To) 


|0| * At 

At ~- , 

n — 1 ATq 


where At = T — To and 0 is the time-average of 0 during At. 
ConsideringEq.(63), thus, it follows that ^ q, so that, 

/(t)-^ /"£>(t") • V dr" 

n — 1 Jto 

~ -AxgD-V 

~ AxqAxD ■ Df . 


Here in the second step of estimations, we have taken into 
account the time-integral is overwhelmed by the contribution 
in the final period of order Atq due to the factor e''^^ along 
with Atq ^ At; in the final step of estimations, Eq.(64) has 
been used. 

Since we are considering the case / < 0, it follows 
0<O<=^ |/(t)| >ATeAT|£)-D/| . 


We have thus found out the following fact: / and D ■ Df are 
much more important than ^ in the later phase of gravita¬ 
tional contractions and that when / is negative and |/(t)| is 
large enough to overwhelm \ D- Df\, the contraction rate gets 
larger. 

Now we turn to Eq.(51), which with the help of Eq.(62) can 
be recast as 


de\ 

dx ) 


= 2 


02 \ df 

-Aa.a /-20^ 


n — 1 


20 


^£)-V-2(7-JDv-(f?..)° 


(67) 


where we have set cOah = 0 by imposing the “twist-less” con¬ 
dition Eq.(49). 

Let us search for the condition for (^) < 0. In later phase 
of the geodesic congruence contractions, the threshold condi- 
tion (^) =0 becomes 


dx 


I0| 


n — 1 
1 


aa 


/ 




It is quite impressive that Eq.(68) is qualitatively similar to 
Eq.(66) as differential equations. Looking at the first two 
equations in Eq.(48), one reasonable estimation is |0| > \(yab\ 
as far as this inequality holds initially. Provided that the 
strong energy condition is preserved, the value of | (/?..)“ | 
should not be very large and it is a reasonable estimation that 
|0| ^ I (/?..)“ I. Thus we reach almost the same conclusion 
that / and D ■ Df are important in the later phase of grav¬ 
itational contractions and that, when / is negative, the large 
|/(t)| (compared to \D- Df\) gives rise to the large contrac¬ 
tion rate. 


VII. SUMMARY AND DISCUSSIONS 


In this paper, we have attempted to construct a mathemat¬ 
ical framework for analyzing the later stages of gravitational 
contractions in terms of the focusing properties of timelike 
geodesic congruences. 

The framework has been formed as a combination of 
a key equation, which relates the geometry variations to 
the geodesic shifts, and a set of equations, which relate 
the changes in the geodesic convergence properties to these 
geodesic shifts. Then, what we have constructed can be 
viewed as a set of equations describing the linear response 
of the geodesic convergence properties to arbitrary geometry 
variations. Since the geometry variations can be arbitrary and 
their origins need not be specified, we might be able to use 
them as probes to investigate the gravitational contraction pro¬ 
cesses in the present framework. 

It has then turned out that the equations get simplified dras¬ 
tically in the case of the conformal variations. We have then 
studied the latest phase of gravitational contractions in the 
case of conformal geometry variations, and have found out 
that in the final stage, / and D ■ Df are much more impor¬ 
tant than and furthermore that the contraction rate gets 
larger when / is negative and |/| is large enough to overwhelm 
\D-Df\. 

It shall be fruitful to change the classes of geometry vari¬ 
ations and see the differences in the focusing properties of 
geodesics, as has been done to some degree in Sec.V by re¬ 
stricting the class to the one of the conformal variations and 
in Sec. VI by further restricting the class to the one satisfying 
the strong energy condition. 

The fact that the case of the conformal variations has re¬ 
sulted in enormous simplifications in equations might imply 
that the non-conformal components of geometry variations 
contain important information on gravitational contractions 
which require further investigations. On the other hand, how¬ 
ever, it means that we can make full use of the simplifications 
for the case of the conformal variations, which itself is impor¬ 
tant and intriguing. 

One possible application related to the conformal geometry 
variations might be the application to the black-hole space- 
times, where the event horizon can be characterized as the 
set of the zero-points of the conformal factor in the metric 
by choosing an appropriate system of coordinates. Then the 
conformal geometry variations in this case correspond to the 
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various deformations of the event horizon so that investigating 
their influence on the convergence properties of the geodesic 
congruence might be of significance. 

As a more mathematically oriented application of the con¬ 
formal variations, there are issues regarding the properties of 
the conformal mappings. Studying mathematical properties of 
the Penrose diagrams are a typical example in this category. 

It is expected that the present framework serves as a tool to 
understand more about the late stages of gravitational contrac¬ 
tions. 
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Appendix A: Useful formulas 


Let Yab be a symmetric covariant tensor. Then it follows 

yab = 7^^-2^a7-^+7..^a^b , 

= lab~'^^^yb-y-^‘^^b , (A5) 

where 7 * is the shorthand notation for the spatial part of 7 ^,. = 
7 bc^‘^- We also note the difference between jb 7b the 
sign difference in the first and the second lines. 

Taking the derivative on both sides of Eq.(A5) with 
some straightforward computations, one can obtain another 
formula for jab 7 

^b7cd = ^b7 + 2 ^- (7 B)-^ b-^L^db- 2 -^- ^bTj 

-2{Y_B)b^c^d + ^b7-^c^d + '2-7-^-B-3b > (A6) 

which can most easily be derived by applying Eq.(A3) to UaUb- 
In particular, for a symmetric spatial tensor 7 ^^, Eq.(A 6 ) 
simplifies to 


Below we enumerate useful formulas frequently used in the 
body of the paper. 

Some among them are formulas regarding the ((n — 1) -f 1 )- 
decomposition of several quantities. Eor their derivations, it is 
helpful to use the relation Eq.(2) and its variations, such as 

= . 


^b7_,, = ^b7_,, + ^^c{7_Bhb ■ (A7) 

There are some useful formulas for the second derivatives 
of a function /; 

^Daf = Da‘^-B^aDcf , (A 8 ) 


With the help of these relations, any vector Ua can be de¬ 
composed as 

Ua = U^-U-^a , (Al) 

where u- = Ub^^ (see Sec.IIC for notations employed in this 
paper). 

Eor a scalar function /, then, the application of Eq.(Al) to 
a vector V^/ yields 


and 

^a^bf = DaDbf-Bba^-2^~D^f + Ub^ , 

(A9) 

and also 

DaDbf = h^ah‘^b^c^df + Bba^ ■ (AlO) 


Vaf^Daf--f^a , (A2) 

dz 

where Da is the spatial derivative operator induced from 
(Eq.(7)). 

Eor notational brevity, it is convenient to introduce 


V, = £> 


a 



It is important to note that coincides with the standard V^, 
only when the former is applied to functions (Eq.(A2)). 

Eor a general spacetime vector u‘‘, it follows 


ybu"^ = ^bu‘‘ + ^‘‘{uB)b-^‘‘^bU--u-B\ ■ (A3) 


Eor a spatial vector then, Eq.(A3) reduces to 

Vbu‘‘ = ^bu‘‘ + ^%uB)b . (A4) 


It is obvious that the contravariant index a in Eq.(A3) and 
Eq.(A4) can freely be lowered and the corresponding formu¬ 
las with the covariant index a also follow. 


By the contraction of Eq.(A9), we get 

df d^f 

Af = D.Df-eJ^-^ . (All) 

On the other hand, by the anti-symmetrization of Eq.(A9), we 
get 


{DaDb — DbDa)f = —2(Oab— , (A 12 ) 

dz 

implying that the induced covariant derivative Da is torsion- 
free iff (Oab = 0. Thus the following statements are mutually 
equivalent: 

( 1 ) (Oab = 0 ^a^b = ^b^a DaDbf = DbDaf- 

(2) The timelike geodesic congruence in question is hy¬ 

persurface orthogonal, i.e. can be foliated by smooth 
(n — 1 )-dimensional orthogonal sections. 

(3) Bab coincides with the extrinsic curvature Kab for the 

smooth orthogonal sections of 
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We here pay attention to a useful formula including the o- 
derivative. By taking the o-derivative on both sides of h^t = 
gab + or /i"'’ = g"'’ + and using Eq.(12), Eq.(13), 
Eq.(22), Eq.(Al) and Eq.(A5), we get 

h ab 'Y_ab 5 

^= -Jab + (Vfo + 7j) • (A 13) 


connection 1 -forms as 

A\ = A'^o = l&‘‘ (fe= 1,2,3) , 

A'™ = -A“i=- 2 ^ 0 ” (m=l, 2 ) , 
art 

A^^ = -A^2 = -—coi^®^ ■ 

ar 


(B4) 


Taking the o-derivative on both sides of hj’ = h/hc^ and us- From Eq.(B4), we can derive the expression for the Riemann 
ing Eq.(A13), it is straightforward to show that curvature as 


h 


b 

a 


= 0 . 


(A14) 


Appendix B: Basic facts on the Friedmann-Robertson-Walker 
model 

We here summarize the basic points on the Eriedmann- 
Robertson-Walker (ERW) model discussed in Sec.VB. 

The ERW spacetime is a spacetime foliated by maximally 
symmetric (i.e. homogeneous and isotropic) spatial hypersur¬ 
faces. The metric for the ERW spacetime is then given as 

ds^ = —dt^ + a^{t) (r)d+ r^d'd^ + sin^ ’ddcp^) , 

F{r):=^=^ {k =-1,0,1) . (Bl) 

The world-lines of observers “standing still” relative to the 
system of coordinates, Eq.(Bl), is described by x"(f) = 

ro t?o (Po) with ro, tSo % being some constants. They 
are easily shown to be timelike geodesics. Let us call these 
geodesics standard geodesics compatible with the system of 
coordinates, Eq.(Bl). 

The spatial part of the coordinate system ({r, 0,0}), thus, 
is interpreted as the collection of the “names” of the observers 
free-falling along each standard geodesic. We also note that 
the time-coordinate t coincides with the proper-time T for the 
free-falling observer along a standard geodesic. Hereafter 
each standard geodesic is assumed to be parametrized by the 
proper-time T, coinciding with t. Based on these facts, we 
see that the coordinate system {t,r,'d-,(p} forms a co-moving 
coordinate system. 

It is sometimes convenient to shift from the coordinates 
{t,r,'d-,(p} to the local pseudo-ortho-normal frame of coor¬ 
dinates { 0 *^, 0 ^, 0 ^, 0 ^} defined as 

= dt , &^ = a{t)F{r)dr , @^ = a{t)rd'& , 

&^ = a{t)rsm'dd(p , (B2) 

which makes the expression for ds^ simply 

c/s 2 = _(@O)2_^(0l)2_^(@2)2_^(@3)2 ^ ^ 33 ^ 

Let be the tangent vector of a standard geodesic and let 
hab be given by Eq.(2). We then see that <^“ = ^(1 00 0), 
1 0 0 0 ) and h^b = diag (0 1 1 1 ) in the { 0 "} frame. 

Following the standard procedures for differential 
forms [15, 16], then, we first get the expression for the 


= - (^ = 1,2,3) , 

2 pi 


^ mlm — 1 

a 




^ 323 = - - 


dV 


1J a^r^F^ ’ 


(B5) 


with the other independent components being zero. From 
Eq.(B5), we get useful coordinate-independent relations 


a j. a 

^abc- — '^~^achb , ^a-b- — ^ab 
a - - a 


Then the Ricci curvature becomes 
, d 


(B 6 ) 


Rqo = — 3 - 


Rn=2 


a IF' 


Rll =^33 = 2 - -f - 


a o^rF^ ’ 

^ a rF'+F^-F 


a 


a^j-2pT> 


(B7) 


with the other components being zero. From Eq.(B7), we get 
other useful coordinate-independent relations 


Ra. = 2>~^a , R.. = 3— . 
a a 


Finally the scalar curvature becomes 


R = 6--f6 
a 


+ 2 


2rF'+F^-F 

a2j.2p3 


(B 8 ) 


(B9) 


For reference, let us write down the expression for the Einstein 
tensor also; 


Goo = 3 


Gii=- - -2- + 


2rF'+F^-F 

a^j-2p3 

a 1 -F^ 


g2j2p2 


G 22 = G 33 = — 


- 2 - 


with the other components being zero. 


(BIO) 
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